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Abstract. Spin-dependent transport through a two-level quantum dot in the 
sequential tunneling regime is analyzed theoretically by means of a real-time 
diagrammatic technique. It is shown that the current, tunnel magnetoresistance, 
and shot noise (Fano factor) strongly depend on the transport regime, providing 
a detailed information on the electronic structure of quantum dots and their 
coupling to external leads. When the dot is asymmetrically coupled to the leads, 
a negative differential conductance may occur in certain bias regions, which is 
associated with a super-Poissonian shot noise. In the case of a quantum dot 
coupled to one half-metallic and one nonmagnetic lead, one finds characteristic 
Pauli spin blockade effects. Transport may be also suppressed when the dot levels 
are coupled to the leads with different coupling strengths. The influence of an 
external magnetic field on transport properties is also discussed. 



PACS numbers: 72.25.Mk, 73.63.Kv, 85.75.-d, 73.23.Hk 
1. Introduction 

Quantum dots have already paved their way to become underlying devices of 
magnetoelectronics and spintronics [U O [3l [4] - not only because of beautiful physics 
emerging in those systems, but, more importantly, due to possible future applications 
and due to the possibility of manipulation of a single spin [SI [HI El [5]. Transport 
characteristics of quantum dots coupled to nonmagnetic leads have already been 
extensively studied both theoretically and experimentally. In nonmagnetic single- 
electron devices it is mainly the electron charge which determines the system transport 
properties [9l [TOl [HI [12] . If, however, a quantum dot is coupled to ferromagnetic leads, 
transport characteristics strongly depend on the spin degree of freedom, leading for 
example to the suppression of current when the alignment of magnetic moments of 
the leads switches from parallel to antiparallel [13l [HI [15] . The difference between the 
currents flowing through the system in these two magnetic configurations defines the 
so-called tunnel magnetoresistance (TMR). 

Theoretical considerations of electronic transport through a quantum dot weakly 
coupled to ferromagnetic leads were restricted mainly to single-level quantum dots. 
Transport properties of such systems were analyzed in the sequential tunneling regime 
[HI [HI [17] , as well as in the cotunnehng regime [HI [HI [SO]- In real systems, however, 
usually more than one energy level participate in transport, leading to more complex 



Transport through two-level quantum dots weakly coupled to ferromagnetic leads 2 

and interesting transport characteristics [HJ [551 123] • Several experimental realizations 
of quantum dots attached to ferromagnetic contacts have already been reported, 
which include self-assembled dots in ferromagnetic semiconductors [24j . ultrasmall 
metallic grains [3 , granular structures [55] , carbon nanotubes [571 UHl HE ISO] , 
single molecules [31], or magnetic tunnel junctions [35] , 

In this paper we consider transport through two-level quantum dots coupled to 
ferromagnetic leads, and restrict our considerations to the sequential tunneling regime. 
The dot is described by the Anderson-like impurity Hamiltonian. Thus, the analysis 
may also hold for some magnetic impurities and molecules. With the aid of the real- 
time diagrammatic technique, we calculate the current /, differential conductance G, 
and Fano factor F in the parallel and antiparallel magnetic configurations, as well as 
the corresponding TMR. The Fano factor, F = S/Sp, describes the deviation of the 
zero- frequency shot noise S from the Poissonian shot noise Sp — 2e\I\ (corresponding 
to uncorrelated electronic transport). 

In the following we analyze transport properties in the two situations: (i) when 
the dot is symmetrically coupled to the leads, and (it) when the dot is coupled 
asymmetrically to the leads. In both cases transport characteristics are shown to be 
strongly dependent on the transport regime. As a result, variation of TMR with the 
transport and gate voltages is significantly different from that for single-level quantum 
dots. Furthermore, the Fano factor is found to exhibit a nontrivial dependence on the 
magnetic configuration of the system and spin polarization of the leads. For symmetric 
coupling to ferromagnetic leads, the Fano factor becomes divergent in the parallel 
configuration when the leads' polarization tends to unity, while in the antiparallel 
configuration shot noise is Poissonian. On the other hand, in the case of asymmetric 
coupling of the dot to external leads, we find transport regions where current is (partly) 
suppressed. These transport regions are accompanied by NDC and a super-Poissonian 
shot noise. If the dot is coupled to one half-metallic and one nonmagnetic lead, the 
blockade regions are associated with particular occupation of a dot spin state and can 
be referred to as the Pauli spin blockade regions. In addition, we also discuss the effect 
of an external magnetic field on the transport characteristics. The TMR displays then 
a rather complex behavior with the bias voltage. 

We note that such effects as spin blockade, NDC, sub- and super-Poissonian shot 
noise, and tunnel magnetoresistance occur also in single-level quantum dots [151 US] ■ 
Although the physical mechanisms responsible for some of the above effects in two- 
level quantum dots are similar to those in single-level dots, the resulting transport 
characteristics in two-level dots are much more complex. The analysis of these effects 
in the case of two-level quantum dots coupled to ferromagnetic leads is the main 
objective of the present work. 

The paper is organized as follows. In section 2 we describe the model. Method 
used to calculate transport characteristics is presented briefly in section 3. Results of 
numerical calculations are shown and discussed in section 4, where we first consider 
the case of a quantum dot coupled symmetrically to the leads, and then analyze 
transport through quantum dots with unequal couplings. Section 4 also includes 
some approximate analytical expressions obtained for the most characteristic transport 
regions. The corresponding formulas may be useful in discussing future experiments. 
Finally, summary and conclusions are given in section 5. 
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2. Model 



The system considered in this paper is shown schematically in Fig. [T] and consists of 
a quantum dot with two orbital levels coupled to ferromagnetic leads. The net spin 
moments of the leads are assumed to be collinear, i.e., they can form cither parallel or 
antiparallel magnetic configuration, as indicated in Fig. [TJ The Hamiltonian H of the 
system includes four terms, H = Hi^ + Hj^ + Hu + Ht- The first two terms describe 
noninteracting itinerant electrons in the leads, = X^ko- ^rkacli^^Crka for the left 
(r = L) and right (r = R) leads, with e,ko- being the energy of an electron with the 
wave vector k and spin cr in the lead r, and cl^^^ (crka) denoting the respective creation 
(annihilation) operator. The quantum dot is described by the following Anderson-like 
Hamiltonian: 

ja j aa' j 

where dj^ {dja) is the creation (annihilation) operator of an electron with spin a on 
the jth level {j — 1,2), Sj is the corresponding single-particle energy, and Uj^r is the 
particle number operator, Ujcr = d^-^dja- The on-level Coulomb repulsion between 
two electrons of opposite spins is described by JTj, whereas the inter-level repulsion 
energy is denoted by U' . The forth term in Eq. ([1]) describes the Zeeman energy, 
with A = g^sB being the Zeeman splitting of the energy levels {B is an external 
magnetic field along the magnetic moment of the left electrode). To describe the 
energy structure of the dot we introduce the level spacing Se ~ 62 ~ £1 and define 
El = £. Since the different orbital levels can couple differently to an applied magnetic 
field, 5e can be controlled by changing the magnetic field [33], while e can be tuned 
by a gate voltage. 

Tunneling p ribed by the 




Figure 1. (color online) Energy diagram of a two-level quantum dot coupled 
to ferromagnetic leads. For clarity reasons the energy diagram is shown here 
for A = (7' = 0. However, in the paper we consider the case when the Coulomb 
interaction between electrons occupying different orbitals is described by a nonzero 
value of the parameter U' . The leads' magnetizations can form either parallel or 
antiparallel configurations. The arrows indicate the net spin of the leads. 
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Table 1. The dot eigenstates \x) and their energies e^- 
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Hamiltonian, 

^T= ^ I](^r,4k.rfj.+i*,4-C'-k.) , (2) 
r— L,R kjcr 

where trj denotes the tunnel matrix elements between the lead r and the jth dot level. 
Couphng of the jth level to external leads can be described by V^j = 2T:\tj.j\^ , with 
jO^ being the spin-dependent density of states in the lead r. By introducing the spin 
polarization of the lead r, = (p^ — /9~)/(p+ the coupling parameters V^^ can 

be expressed as T^- ^ = Tr]{l ± Pr), with Trj = (r+ + T^^)/2. Here, T^^ and T^^ 
describe the coupling of the jth level to the spin-majority and spin-minority electron 
bands, respectively. In general, each dot level may be coupled to the leads with a 
different strength. Moreover, the coupling strengths may be energy dependent. In 
this work, however, they are assumed to be constant within the electron bands. For 
example, in the case of experiments performed in the weak coupling regime by Kogan 
et al [31] , the coupling strength was found to be of the order of tens of /ueV. 

3. Method 

We analyze the spin-polarized sequential transport through a two-level quantum dot. 
The first-order tunneling gives the dominant contribution to charge current for voltages 
above a certain threshold voltage and is exponentially suppressed in the Coulomb 
blockade regime. The effects due to higher-order tunneling, e.g., cotunneling p!8l |35] 
are not included. Taking into account the two orbital levels of the dot results in sixteen 
different dot states |x) (|x) = Ixi)lx2), with \xj) corresponding to the j-th level), as 
listed in Table [1] The system is symmetrically biased and we assume equal capacitive 
couplings to the left and right lead, so the dependence of the dot energy levels on the 
bias voltage may be neglected. The energies of the corresponding dot eigenstates are 
listed in Tabled 
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Transport is calculated with the aid of the real-time diagrammatic technique 
[551 [571 [T51 HT]. The technique consists in a systematic perturbation expansion of 
the dot density matrix and the current operator in the dot-lead couplings Trj. In 
order to calculate the stationary occupation probabilities, current, and shot noise in 
the lowest order in the dot-lead coupling, we follow the procedure developed in Refs. 
[STlET], and introduce the respective self-energy matrices: W, W'-, W''. The ma- 
trix W contains self-energies with one arbitrary row xo replaced by (F, . . . , F), which 
is due to the normalization of the probabilities, J^x^x — 1- The elements W^x' of 
the matrix W describe the first-order tunneling transitions between the many-body 
dot states |x) and \x')- They are given by [21], Wxx' = , + W^ ,, where , = 



for 



X + x', while = -Ex'^x^x'x' ^ith /.(e) = l/[e(^-^^r)/k^T ^ ^ 

being the electrochemical potential of lead r. The second matrix, W'-, denotes 
the full self-energy matrix with one external vertex, resulting from the expansion 
of the tunneling Hamiltonian, replaced by the current operator. Finally, the third 
matrix, W'', consists of self-energies with two external vertices replaced by the 
current operator. The current operator / is defined as / = (/r — with 

= Ekcr Ej {^rjC^r^icr'^o'y - K^'^ja^rV.c^ being the current flowing from the 

lead r to the dot. The elements of the matrices W' and W'-'- can be expressed 
in terms of l^^x' as [21], W\^, = [e(iVx' " ^x) " ©(^x " ^x')] (^^xx' ~ ^xx') 
^xx' ^ ^ 2(5^x')^xx'/4, respectively, where Nx = J^ja^jcr and 6 (a;) is the step 
function. 

In this paper we calculate transport in the first-order with respect to tunneling 
processes. In the collinear configurations discussed here, one might also apply a simpler 
method based on the master equations. However, we apply a more sophisticated 
technique, which enables easy extensions to noncoUinear magnetic configurations and 
facilitates including the influence of second- and higher-order tunneling processes 
|38| . Furthermore, the real-time diagrammatic technique employed here allows one 
to calculate the current and shot noise in a fully systematic way order by order in the 
dot-lead coupling strength. 

Having calculated the respective matrices, the stationary probabilities can be 
determined from the following equation, 

(WP)^ = TSx^, , (3) 
where P is the vector containing the occupation probabilities. In turn, the sequential 
current flowing through the system can be calculated from 

I=^Tr{W'P}, (4) 

with Tr{A} denoting the trace of the matrix A. On the other hand, the zero-frequency 
shot noise, S ^2 dt ((/(<)/(0) + /(0)/(i)> - 2(/)2^ , is given by 

2 

S = — Tr{W"P -I- W'PW^P} , (5) 

h 

where P is determined from the equation WP ~ Q, with Qx'x ~ 

(P^,-Sx'x)il-Sx'xo)- 

Knowing the current / and the zero-frequency shot noise S, one can determine 
the Fano factor F, F = 5/(261/1). The Fano factor describes the deviation of S from 
the Poissonian shot noise given by Sp = 2e\I\. 
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4. Results and discussion 

In this section we present numerical and analytical results on the current /, TMR, 
and Fano factor F for the two-level quantum dot weakly coupled to the leads. The 
TMR effect is phenomenologically defined as [131 ITS] 

TMR = ~ , (6) 
Iap 

where Ip (Iap) is the current flowing through the system when magnetic moments of 
the leads are aligned in parallel (antiparallel) . 

In our considerations we will distinguish between quantum dots coupled 
symmetrically and asymmetrically to the leads. The asymmetry may be introduced 
in various ways - for example by attaching the quantum dot to leads with different 
spin polarizations (e.g., to one half- metallic and one nonmagnetic lead). On the other 
hand, if the leads have the same spin polarizations, their coupling to the dot levels may 
be different. In the following we first consider the case of quantum dots symmetrically 
coupled to the leads and then we shall proceed to analyze transport through quantum 
dots asymmetrically coupled to the leads. In order to simplify the following discussion 
of numerical and analytical results, we set Ui — U2 — U' = U . 



Quantum dots symmetrically coupled to ferromagnetic leads 

In Fig. [2] we show the differential conductance as a function of the bias voltage V 
and level position e. Experimentally, the position of the dot level can be tuned by a 
gate voltage, therefore Fig. [2] effectively shows the bias and gate voltage dependence 
of the differential conductance. The upper part [Fig. [U^a)] corresponds to the parallel 

the e, 

whi :e 

in £ ;o 

the le 



G (e^/h) 






eV/U 



Figure 2. (color online) The differential conductance G = dl/dV as a function 
of the bias voltage and level position in the parallel (a) and antiparallel (b) 
magnetic configurations for the parameters: k^T = F, 5e = 25r, U = SOF, 
A = 0, PL = PR P = 0.7, and F,.j = F/2 (r = L, R, j = 1, 2). 
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Figure 3. (color online) The tunnel magnetoresistance as a function of the bias 
voltage and level position for the same parameters as in Fig. [2] 




conductance when the leads' magnetizations are antiparaUel. The majority (minority) 
electrons of the source electrode tunnel then to the minority (majority) electron band 
of the drain electrode, leading effectively to a conductance smaller than that in the 
parallel configuration. If the spin asymmetry in the density of states increases (the 
spin polarization of the leads increases), the difference in conductance between these 
two magnetic configurations becomes larger, which leads to an increase in TMR. 

Another feature visible in Fig. [2] is the diamond-like structure of the differential 
conductance, which is associated with discrete charging of the quantum dot. The 
diamonds around = correspond to the Coulomb blockade regions. When lowering 
position of the dot levels, the charge of the dot changes successively. More precisely, 
the dot is empty for e >0, occupied by one electron for > e > —U, doubly occupied 
for -[/>£> -(2C/+(5e), occupied by three electrons for ~{2U+Se) >e> -(3C/+fc), 
and the two orbital levels of the dot are fully occupied for —{3U + Se) > e. In all 
these transport regions the dot is in a well-defined charge state, and the sequential 
tunneling is exponentially suppressed [3 H]. If the bias voltage is increased above 
a certain threshold voltage, the current starts to flow due to first-order tunneling 
processes. When the charging energy of the dot is much larger than the thermal 
energy, one observes then a well-resolved step in the current as a function of the bias 
voltage. In the density plots shown in Fig. [2l this can be seen in the form of lines that 
clearly separate the Coulomb blockade regions from transport regions associated with 
consecutive charge states taking part in transport. When the bias voltage increases 
further, additional steps (and consequently lines in Fig. [5]) arise at voltages where new 
states becomes active in transport. 

Figure [3] presents the corresponding TMR as a function of the bias and gate 
voltages. By comparison with Fig. [2l one can easily identify the regions of Coulomb 
blockade. First of all, it is worth noting that, depending on the transport region, 
TMR takes several well-defined values ranging roughly from 0.2 to 1.6, see Fig. [3l 
Such behavior of TMR is significantly different from that for a single-level quantum 
dot, where TMR in the sequential tunneling regime acquires only two values [18j . As 
follows from Fig. [3l TMR in the linear response regime is independent of the gate 
voltage. More precisely, it is given by p^/(l — P^), which is exactly equal to a half 
of the TMR calculated within the Julliere model [13]. We recall that the JuUiere 
value of TMR, TMRo = 2p'^ / (1 — p^), is characteristic of a single-barrier planar tunnel 
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junction. In the case of double-barrier planar junctions with ferromagnetic leads, the 
TMR in the sequential transport regime is equal to half of the Julliere value. This 
applies also to sequential transport through quantum dots coupled to ferromagnetic 
leads. However, one should bear in mind that transport in the Coulomb blockade 
regime is dominated by higher-order tunneling, which can considerably modify the 
corresponding TMR. For instance, it was shown recently for single-level quantum dots 
that TMR in the cotunneling regime exhibits a strong dependence on the number of 
electrons on the dot, and for empty or fully occupied dots reaches the JuUiere's value 

Figure [3] also shows that, when increasing the bias voltage V and keeping 
constant position of the dot levels, TMR acquires some specific and well-defined 
values in different transport regions. To analyze this behavior in more details, we 
show in Figs. SJa) and EJa) the bias voltage dependence of the current and TMR 
for the case when the dot level is above {e = U/2) and below (e = —7U/5) 
the Fermi level of the leads at equilibrium. As one can see, the current in both 
magnetic configurations and the associated TMR exhibit characteristic plateaus which 
correspond to different transport regions. When assuming the zero-temperature limit, 
one can derive approximate formulas describing transport in each region. However, 
the analytical solution in a general case is rather complicated due to sixteen different 
eigenstates involved. Therefore, in the following we only present the analytical results 
for the transport regions, where the corresponding formulas are relatively simple and 
transparent. These transport regions are marked in Figs.SJa) and Elja) by consecutive 
numbers. An exemplary analytical calculation of the current in parallel configuration 
for plateau (3) is described in the appendix. 

The analytical expressions for the current /p(AP) and TMR are listed in Tabled] 
First of all, the current depends on the spin polarization of the leads effectively only in 
the antiparallel configuration, while in the parallel configuration it is independent of p. 
Such a behavior is a consequence of the left-right symmetry for each spin channel in the 
parallel configuration, and equal occupation of the charge states active in transport. 
This leads to the independence of the average charge and spin of the dot on the 
polarization factor p, which in turn leads to the charge current independent of p. 
This is not the case in the antiparallel configuration, where the symmetry is absent 
and the current does depend on the spin polarization of the leads. The associated 
TMR disappears for p = and diverges as p — > 1 (this is because Iap for 
p — > I). For a finite spin polarization p, one finds a considerable enhancement of 
TMR in the regions (2), (3), and (4) as compared to the other transport regions. 
The maximum TMR occurs in the region (3) and is given by 4/5 x TMRq. This 
enhancement can be accounted for by realizing that in this transport region all the 
one-particle states take part in transport, leading to an enhancement of the spin 
accumulation in the antiparallel configuration. The doubly occupied states do not 
take part in transport due to the Coulomb energy. In the antiparallel configuration, 
spin-down electrons of the source electrode can tunnel relatively easily to the dot, 
where, however, they have to spend a longer time before tunneling further to the 
drain electrode. Consequently, the current is mainly determined by the tunneling 
probability through the barrier between the dot and the drain electrode. This means 
that a new channel for tunneling, which becomes active when going from the region (2) 
to the region (3), has no significant influence on the current, as it is clearly visible in 
Fig. mja) . In turn, in the parallel configuration this new transport channel contributes 
to the current, increasing this way the difference between the parallel and antiparallel 
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Figure 4. (color online) The current (a) in units of Iq = eF/h and Fano factor 
(b) in the parallel (solid line) and antiparallel (dashed line) configurations as well 
as tunnel magnetoresistance (a) as a function of the bias voltage for e = U/2. 
The other parameters are the same as in Fig. [2] 

configurations (and consequently the associated TMR). 

Although the current flowing through the system in the parallel configuration is 
independent of the spin polarization in transport regions marked by the numbers, we 
note that Ip slightly depends on p in the other regions. This is mainly due to the 
fact that in these regions the occupations of the many-body dot states being active 
in transport are not equal, leading to a small spin accumulation even in the parallel 
configuration. For example, in the region (2) there are 3 states taking part in transport, 
with the occupation probabilities equal to 1/3; in the region (4) there are 11 states 
(one empty, four single-particle and six two-particle states), all equally occupied with 
the probability 1/11, etc. In these regions, the charge current does not depend on 
spin polarization. However, there are also regions where the occupation probabilities 
depend on the spin polarization of the leads, such as for example the region between 

(3) and (4) . In this region 10 dot states participate in transport [the same as in region 

(4) except for the state with doubly occupied second level] . We find that in this region 
the dot states do not take part in transport on an equal footing, which leads to small 
spin accumulation. This is a consequence of the interplay between the energy structure 
of the dot and ferromagnetism of the electrodes. The spin accumulation disappears 
when transport goes through single-particle states. This simply follows from the same 
spin asymmetry of tunneling through the left and right barriers. If we have transport 
through doubly occupied states with Coulomb interaction, the situation becomes more 
complex. The Coulomb interaction in some situations may suppress this symmetry, 
leading to spin accumulation. We note that similar spin accumulation for parallel 
magnetic configuration in symmetric junctions was also observed in different systems, 
see eg. [39] , 
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Figure 5. (color online) The bias dependence of the current (a), the Fano factor 
(b) in both magnetic configurations and the TMR (a) calculated for s = —7U/5. 
The other parameters are the same as in Fig. |4] Current is plotted in units of 
lo = eV/h. 

A new and interesting behavior of TMR may be observed in the blockade region, 
where the dot is doubly occupied, see Figs.[2]and[5fa). Contrary to the case of a single- 
level quantum dot, where sequential TMR is constant in the whole blockade regime, 
the TMR displays now a strong dependence on the bias voltage. Once the voltage 
is increased to around eV « U/2, TMR drops rapidly from 1/2 x TMRq to about 
1/10 X TMRq and then, when further increasing the bias voltage, TMR increases to 
the value corresponding to the plateau (7), see Tab.^l This leads to a deep minimum 
in TMR clearly visible in Fig.[5Ua), which can be also seen in Fig. [3] as the two black 
areas in the middle of the figure. 

We have also calculated the Fano factor Fp and Fap in the parallel and 
antiparallel magnetic configurations, respectively. The zero-frequency shot noise can 
be then found from the knowledge of the current / and the Fano factor F. The 
Fano factor for the two situations discussed above, and corresponding to Figs. IH^a) 
and[5l^a), is shown in Figs.|4Ub) and[5fb) for both magnetic configurations. Similarly 
to the case of current and TMR, the Fano factor acquires roughly constant values, 
different in different transport regions. The relevant analytical formulas, obtained with 
the same assumptions as before, are given in Table O Due to the spin asymmetry in 
the coupling of the dot to external leads, the bias dependence of the Fano factor 
is significantly different from that in the corresponding nonmagnetic situations |21j . 
For both magnetic configurations of the system, the Fano factor depends on the 
polarization factor p (differently in the two configurations, in general), except for 
the case of empty dot where the exponentially suppressed transport gives rise to 
Poissonian Fano factor Fp — Fap — 1. In addition, if \eV\ ^ fce^T, the Fano factor 
becomes divergent due to the thermal noise, which dominates in this transport regime; 
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Table 2. Analytical expressions approximating the current 7p (/ap)i the Fano 
factor Fp (-Fap) in the parallel (antiparallel) magnetic configuration, and the 
TMR in different transport regimes corresponding to plateaus shown in Figs. |4] 
and |5]for pL = PR = p- The current is expressed in the units of Iq = eV/h. 
In the Coulomb blockade regime, where the current is suppressed, the TMR was 
determined from the linear conductance. 

Plateau /p /ap TMR Fp Fap 
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in the case of y = 0, the noise is given by S — 4fcBTG''", with G''" being the 
Hnear conductance, leading to a divergency of the Fano factor (finite S for 7 = 0) 
[401 I41j . Furthermore, in some transport regions we find Fp > Fap, while in the 
other ones Fp < Fap . As one can see from the expressions listed in Table [21 the 
ratio Fp / Fap depends on the spin polarization of the leads p. For example, in the 
region (2) Fp > Fap only if p > 0.37. Furthermore, if the leads are half- metallic, 
the Fano factor in the parallel configuration diverges as p ^ 1, except for regions 
(1) and (6), where Fp = 1 and Fp — 1/2, respectively. On the other hand, in the 
antiparallel configuration the Fano factor tends to unity for p ^ 1, except for the 
Coulomb blockade regime with two electrons trapped in the dot, where we find super- 
Poissonian shot noise, Fap « 3, see Fig. [Sjjb). 

The super-Poissonian shot noise in the Coulomb blockade regime with two 
electrons trapped in the dot, shown in Fig. [5ljb), can be accounted for as follows. 
In the ground state the dot is occupied by two electrons on the first level (of lower 
energy) . Assume the voltage corresponding to the maximum of the peak in the Fano 
factor, eV/U ~ 0.5. The system is still in the blockade regime [compare Fig. [SJa)], 
where the current is exponentially suppressed. A nonzero small current can flow due 
to thermal excitations. There is an exponentially small probability that one electron 
leaves the dot and the other one jumps to the dot either to the same (first) or to the 
second energy level (of higher energy) . If it tunnels to the second energy level, then 
transport through this level is allowed and electron can easily leave the dot, while 
another one can jump to the same level or to the level of lower energy. If it tunnels 
to the same level, further tunneling processes are allowed. If it tunnels to the low 



Transport through two-level quantum dots weakly coupled to ferromagnetic leads 12 

energy (first) level, the system is blocked again. This leads to large fluctuations in 
the current, and consequently to super-Poissonian shot noise shown in Fig. [SUb). This 
behavior does not results from magnetism of the electrodes and persists even in the 
case of nonmagnetic leads, where the Fano factor is approximately F 3. The effect 
is only quantitatively modified by ferromagnetism of the electrodes - the magnitude of 
the super-Poissonian shot noise is different in parallel and antiparallel configurations, 
as shown in Fig. [Hb). 

On the other hand, the increase of the Fano factor in the parallel configuration, 
observed when p — > 1, is due to the enhanced spin asymmetry in transport processes 
through the dot [T6l|42]. Consider for instance the region (2) in Fig. [IJb). Spin-up 
electrons tunnel then relatively easily to and out of the dot. However, when a spin- 
down electron appears on the dot, it blocks transport for relatively long time, leading 
to large fluctuations in the current. These fluctuations increase when p ~f 1, giving 
rise to super-Poissonian Fano factor, which is characterized by a divergent component 
Fp ~ (1— In turn, such fluctuations are absent in the antiparallel configuration 
and the shot noise is Poissonian, Fap ~ 1 for p — > 1. 

4-2. Quantum dots asymmetrically coupled to the leads 

In the previous subsections we have analyzed transport through quantum dots coupled 
symmetrically to ferromagnetic leads (equal spin polarizations of the leads and equal 
coupling parameters). In that case transport characteristics were symmetric with 
respect to the bias reversal. However, when the dot is coupled asymmetrically to the 
leads, the I — V curves are no longer symmetric with respect to the bias reversal, 
leading to further interesting effects. In the following, we discuss transport properties 
of quantum dots asymmetrically coupled to the leads and consider two special cases. 
The first case concerns the quantum dot coupled to one nonmagnetic and one half- 
metallic lead (pl = P, Pr = and Fi^ = r2r')) whereas in the second case the 
dot is coupled to the leads with equal spin polarizations but with unequal coupling 
strengths (pl — Pr = P and Fi^ ^ ^2r')- It ha-s been shown recently that transport 
characteristics of a single-level quantum dots coupled to one half-metallic and one 
nonmagnetic lead display a diode-like behavior [Kl UHl 133] • On the other hand, in 
the nonmagnetic limit of the second case a negative differential conductance has been 
found [21]. 

4-2.1. Unequal spin polarizations of the leads When one of the leads is half-metallic 
{p = 1) and the other one is nonmagnetic {p = 0), transport characteristics become 
asymmetric with respect to the bias reversal. Furthermore, the current can be 
suppressed in certain bias regions, and this suppression is accompanied by the 
occurrence of NDC. This basically happens when the electrons residing in the dot 
have spin opposite to that of electrons in the half-metallic drain electrode. In Fig. [6] 
we show the current and Fano factor for a quantum dot coupled to half-metallic (left) 
and nonmagnetic (right) lead as a function of the bias voltage. For the parameters 
assumed, the current is suppressed in certain ranges of positive bias voltage. As follows 
from Fig. [G] there are three blockade regions. To facilitate the further discussion, we 
label these blockade regions with the numbers, as indicated in Fig.[6l^a). On the other 
hand, for negative bias voltage, the current changes monotonically with the transport 
voltage, as one can see in the inset of Fig. [61(a). 

In the blockade region (1), 1 < eV/U < 2, the dot is in the state | |)|0) and 
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Figure 6. The current (a) in units oi Iq = eF / h and Fano factor (b) as a function 
of the bias voltage for pL = P = 0.95, 0.98, pR = 0, while the other parameters 
are the same as in Fig. [2] The inset in part (a) shows the current in the whole 
range of the bias voltage. 

the spin-down electron residing in the dot has no possibility to tunnel further to 
the left lead, which leads to suppression of the current. The blockade region (2), 
2 < eV/U < 3, is associated with the occupation of the states | i)|0) and |0)| J,). The 
current is then prohibited due to the full occupation of the single-particle spin-down 
states. When increasing the bias voltage further, 3 < eV/U < 4, the blockade of the 
current becomes suppressed [see the plateau between regions (2) and (3) in Fig.[Sl[a)], 
which is due to a finite occupation of state | Ti)|0). Although tunneling of spin-down 
electrons is then blocked, the current is still carried by spin-up electrons. In turn, the 
blockade region (3) occurs for 4 < eV/U < 5, where the dot is in the triplet state 
11)11) and tunneling is also suppressed. Thus, the current is blocked when the total 
dot spin Sz is either Sz — —1/2 or 5*2 = —1, i.e., the spin of electrons on the dot 
is opposite to that of electrons in the half-metallic lead. There is no suppression of 
the current for negative voltage. This is because the electrons residing in the dot can 
always tunnel to the nonmagnetic drain electrode. 

The current blockade in the regions (1) to (3) in Fig. ^a,) for positive bias is 
not due to the charging effects as in the Coulomb blockade regime, but due to a 
particular occupation of the dot spin state. Such blockade is frequently referred to 
as the Pauli spin blockade, and has already been found in single-dot and double-dot 
systems [44l |45l |46]. Here, we show that the Pauli spin blockade can occur in single 
two-level quantum dots, provided they are coupled to half-metallic lead(s). It is further 
worth noting that the blockade of the current in Fig. [DJa) is not complete - there is a 




Transport through two-level quantum dots weakly coupled to ferromagnetic leads 14 



small leakage current inside each blockade region, which resuhs from the fact that the 
assumed spin polarization of the half- metallic lead is not exactly equal to unity. When 
spin polarization is increased, the current in the blockade regions decreases (compare 
the curves for p — 0.95 and p = 0.98). 

The spin blockade of charge current leads to the super-Poissonian shot noise, i.e., 
the corresponding Fano factor is larger than unity. On the other hand, the Fano factor 
outside the spin blockade regions is sub-Poissonian (smaller than unity), as shown in 
Fig. [6ljb). The enhancement of the shot noise in the Pauli blockade regions is a 
consequence of large spin asymmetry in the tunneling processes. The occurrence of a 
spin-down electron on the dot prevents further tunneling processes for a longer time, 
while spin-up electrons on the dot can escape much faster, allowing further tunneling 
processes. This gives rise to large current fluctuations, and consequently also to Fano 
factors much larger than unity. 

Assuming the zero temperature limit, one can find approximate expressions for 
the current (in the units of Iq = eT/h) in the spin blockade regions 

/(I) = (7) 



for the first (1), 



/'^> = ^. ,8) 

5 — 

second (2), and 

(3) ^ 2(l-p2)(2-p2) 

for the third (3) blockade region, see Fig. [DJa). The above formulas show that 
the current vanishes when p ^ 1. This is due to the total blockade of charge 
transport by spin-down electrons in the dot. This can be also concluded from the 
respective occupation probabilities, which are given by -P||^)|o) = (1 — p'^) for 

the first blockade regime, ^'|j^)|o) = ^\o)\i) ^ (1 +P)/(5 ^P^) for the second one, and 
~ +-P)^/(11 ~ + P^) for the third blockade region. It is clearly evident 
from these formulas that the corresponding occupations approach unity for p — > 1, 
leading to the full spin blockade. It is also interesting to note that in region (3) one 
observes formation of a pure triplet state. 

With the same approximations as above, one can find the corresponding formulas 
for the Fano factor in the spin blockade regions 

for the region (1), and 

_ 17-t-30p2+p4 

for the region (2). From the above expressions follows that ==5/9 and 

= 17/25 for p = 0, i.e., the shot noise is sub-Poissonian, However, for p —^ 1 
the noise becomes super-Poissonian, F'^^'> — 3 and F^^^) — 51/16. In the third spin 
blockade regime, the analytical formula for the Fano factor F*^^^ is too cumbersome 
to be presented here. It is however worth noting that F^^^ « 0.55 for p — 0, while 
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Figure 7. (color online) The current (a) in units of /q = er/?l and the Fano 
factor (b) in the parallel (solid line) and antiparallel (dashed line) configurations, 
and the resulting TMR (c) as a function of the bias voltage for e = 25r and 
Fli = Fri = 50Fl2 = = F/2. The other parameters are the same as in 

Fig. El 

^(3) = 5 for p 1. In addition, from the analytical expressions for the Fano factor 
one can easily estimate the value of spin polarization p at which the shot noise becomes 
super-Poissonian. This gives 

> 1 for p > and F^^) > 1 for p > 1/^5. 

4-2.2. Unequal coupling of the dot levels to the leads It has been shown recently 
that a negative differential conductance (NDC) can occur in transport characteristics 
of two-level quantum dots, when the two orbital levels are coupled to nonmagnetic 
leads with different coupling strengths [21] ■ In the following, we consider the dot 
coupled to ferromagnetic leads of equal spin polarizations, but with unequal coupling 
strengths. For one bias polarization we find narrow bias regions where NDC occurs, 
followed by the regions where the current is partly suppressed. The suppression of 
current is accompanied by a super-Poissonian shot noise and may lead to negative 
TMR, as shown in Fig. [3 where the current and Fano factor in the parallel and 
antiparallel magnetic configurations as well as the corresponding TMR are plotted 
as a function of the bias voltage. As one can see, NDC and the associated regions 
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where charge current is partly suppressed occur in both magnetic configurations of 
the system. The suppression regions are marked with (1), (2) and (3), see Fig. [Tl^a), 
and occur for 2 < eV/U < 3, 4 < eV/U < 5, and 6 < eV/U < 7, respectively. 

In the blockade region (1), the dot can be occupied by at most one electron. The 
blocking mechanism is associated with an increased occupation of the second orbital 
level and decreased occupation of the first level. Since the coupling of the second 
level to the left lead (drain for electrons at positive bias) is the weakest coupling in 
the system (the others are assumed to be equal) , the average occupation of this level 
is larger than the occupation of the first level. Thus, an electron that has tunneled 
from the right lead to the second level spends a long time in the dot before it tunnels 
further to the right lead and, consequently, blocks transport through the first level. 
In other words, the onset of electron tunneling to the second level of the dot leads 
to the suppression of current (NDC and the associated suppression region). In turn, 
no NDC and current suppression takes place for negative bias voltage, as now an 
electron that has tunneled to the second level can easily tunnel further to the drain 
(right) electrode. From the above discussion follows that the suppression of current 
for positive bias voltage occurs in those transport regions, where the occupation of 
the second level is larger than that of the first dot level, and the rate for tunneling 
between the second level and the drain electrode is smaller than that for the source 
electrode. The same scenario also holds for the next two blockade regions. In the 
region (2) it is the increased occupation of the state |0)| ti) which is responsible for 
the blockade, whereas in the region (3) transport is mainly mediated through the 
states |0)| ti), I T)l Ti), and | i)| ti)- If the second level is completely decoupled from 
the left lead, transport in the region (3) is still mediated through the single-particle 
states of the first orbital level. This is why there is only a weak suppression of the 
current in this blockade region, as compared to the first and second blockade regions, 
where the current can be fully suppressed. 

In the zero temperature limit one can derive some analytical formulas for the 
current, TMR and Fano factor in all the three blockade regions. First of all, 
it is interesting to note that NDC in the regime (1) occurs when x < 1/3 and 
X < {1 + p'^)/{3 + p^) for the parallel and antiparallel configurations, respectively, 
with X = rL2/r. The parameter x describes the asymmetry between the coupling of 
the second dot level to the left lead and the other couplings, Fli = Fri = Fr2 = F/2. 
Thus, for X ^ the second dot level is completely decoupled from the left lead, 
whereas for x — 1/2 one recovers the symmetric case. The current (in the units 
of Iq) in both magnetic configurations is given by /p^' = 2x/(l + 3a;) and /^p = 
2(1 - p'^)x/[l + 3x + (1 + x)p2], whereas the TMR is 



It is worth noting that the current tends to zero for a; — > 0, whereas TMR approaches 
TMRq. The formulas for the Fano factor are too complex to be written here 
explicitly. However, we find that for a; 0, F^^'^ = (3 + p^)/{l - p^) and 

F^p = (3 + lOp^ — p^)/(l + p^)^ for the parallel and antiparallel configurations, 
respectively, which yields supcr-Poissonian shot noise, irrespective of spin polarization 
of the leads. 

It is also possible to derive analytical formula for the current, TMR, and the 
conditions for NDC in the other blockade regions. However, they are not simple 
and would make the discussion somehow obscure. Therefore, we will only consider 
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two limiting situations, i.e., the symmetric case and the case of a; — > 0, where the 
formulas become relatively simple. When the couplings are symmetric, the blockade 
is lifted and the current, TMR and Fano factor in the blockade regimes are given 
by the expressions corresponding to the third, forth and fifth plateaus, as listed in 
Table [2l When the second level is decoupled from the left lead {x ^ 0), current in the 
second blockade regime tends to zero, while the TMR becomes negative and is given 
by TMR(2) = -V/(3 + p2), see Fig.lHb). 

To understand the physical mechanism leading to the negative sign of TMR let us 
assume a small value of x and that the dot is in the blocked state with two electrons 
on the second level. When a spin-up electron tunnels from the dot to the left lead 
(this tunneling probability is higher than for spin-down electrons), the system becomes 
unblocked, and the next blockade takes place when spin-up electron from the right lead 
tunnels to the dot. This tunneling probability is larger in the parallel configuration 
than in the antiparallel (in the former case it involves spin majority electrons while 
in the latter case spin minority ones). Thus, the blockade is more effective in the 
parallel configuration than in the antiparallel one, leading to negative TMR for small 
values of x, as shown in Fig. EKb). For the parameters assumed to calculate Fig. [71 
the TMR changes sign when x < 0.023. For the Fano factor in the limit of a; ^ 

(2) (2) 

we find F^p > Fp for p < 1. However, if p — + 1 the Fano factor is divergent in 
both magnetic configurations of the system. Furthermore, in the limit of x 0, 
the shot noise becomes super- Poissonian, irrespective of magnetic configuration of the 
system and spin polarization p. For p = 0, the Fano factor is minimum and given by 
Fj,^'> = Ffj] w 4.9. 

Contrary to the first two blockade regions, current in the third blockade region 
is finite for x ^ 0. This is due to the fact that although tunneling through the 
second level is suppressed, current can be still mediated by single-particle states of 
the first level. In the zero temperature limit the current in the parallel and antiparallel 
configurations is given by (in the units of Iq) Ip^ = 1/3 and ifp = (1 — p^)/(3 +p^), 
respectively. This yields the tunnel magnetoresistance TMR(^) = V/(3 - 3p2) ^ 
2/3 X TMRq, which is characteristic for the transport regime where single-particle 
states take part in transport, see the formulas for plateau (2) in Table [H The 
Fano factor in the case of negligible coupling in both magnetic configurations 
is F^^^ = (5 + 3p2)/(9-9p2) and i^|p^ = {5-p^){l + 3p^)/{3+p^)^, respectively. This 
implies that F^p < 1 for all p, while Fp^-' becomes larger than unity for p > 1/V3. 

As follows from the above discussion, a distinctively different transport behavior 
of the system can be found in each blockade region, when the coupling of the second 
level to the left lead is negligible. This is especially visible in TMR which in the first 
blockade region approaches the JuUiere's value, in the second blockade changes sign 
and becomes negative, while in the third blockade region is given by 2/3 of the Julliere 
value. Furthermore, in the two first blockade regions we find a full suppression of the 
current for a; ^ 0, accompanied by a super-Poissonian shot noise. It is also worth 
noting, that the origin of the blockade is now different from that discussed previously 
(Pauli spin blockade), although some qualitative features of the blockade are very 
similar. 
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Figure 8. (color online) The current (a) in units of Iq = eV/h and Fano factor (c) 
in the parallel (solid line) and antiparallel (dashed line) configurations as well as 
the TMR (b) as a function of the bias voltage in the presence of external magnetic 
field, A = sr. The other parameters are the same as in Fig. |2] For comparison 
we also show the TMR in the case of A = 0, see the dotted blue line in part (b). 

4-3. Effect of external magnetic field 

So far we have assumed that the dot levels are spin degenerate. However, in a strong 
magnetic field the Zeeman energy can be large enough (of the order of T or larger) 
to influence transport characteristics significantly [34]. The impact of an external 
magnetic field on transport properties is the most visible in the case of symmetric 
couplings, therefore only this situation is discussed in the following. 

Bias dependence of the current and Fano factor in both magnetic configurations 
and the corresponding TMR are shown in Fig. [5] for the symmetric system and in 
the presence of external magnetic field. Because of a finite Zeeman splitting of the 
dot levels, there is now an asymmetry of transport characteristics with respect to the 
bias reversal. Moreover, NDC occurs now for both positive and negative bias voltages, 
irrespective of magnetic configuration of the system, as clearly visible in Fig. [Ha). The 
NDC is however more pronounced in the parallel configuration than in the antiparallel 
one. 
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The TMR is displayed in Fig. ^h), where, for comparison, the TMR in the 
absence of magnetic field is also shown. If the energy levels are split due to the Zeeman 
energy, TMR exhibits characteristic peaks which occur at voltages corresponding to 
the Coulomb steps in the TV curves. These peaks are more visible for negative bias 
voltage. In addition, a high plateau evolves in the low bias range due to the splitting of 
the dot levels. To determine the analytic formula for the TMR in the linear response, 
we consider the linear conductance, which in the case of A 3> fceT and low temperature 
can be approximated by Gp" ^ {1 + p)/2 and Gap ^ (1 ~ -P^)/2i for the parallel and 
antiparallel configurations, respectively. This yields for the TMR 

TMR"" (13) 

It is worth noting that TMR is now much enhanced as compared to the case of no 
external magnetic field, where the TMR is given by p^/ (1 — p^). Moreover, the linear 
TMR in the presence of magnetic field is even larger than the Julliere value of TMR. 
The enhancement of TMR is due to the fact that when A 3> fceT, the thermally- 
activated sequential transport takes place only through the state | t)|0) and not 
through states | t)|0) and | i)|0) as in the case of A = 0. This considerably increases 
the spin asymmetry of the total current in both magnetic configurations, and gives rise 
to enhanced TMR in the linear response regime. The peak structure of TMR results 
from the corresponding peaks in current, shown in Fig. [Sja). More specifically, when 
the dot levels are spin-split, the two spin contributions to current differ not only in 
magnitude, but also in positions of the steps (peaks in the corresponding differential 
conductance). Adding the two contributions leads to characteristics of the current 
with a pronounced peak structure. This structure is much more visible in the parallel 
configuration than in the antiparallel one. This simply follows from spin asymmetry 
- note that for vanishing Zeeman splitting the two contributions to current become 
equivalent in the antiparallel configuration, but are still significantly different in the 
parallel one. 

The Fano factor in both magnetic configurations is shown in Fig.[5I^c). Contrary 
to the current and TMR, behavior of the Fano factor is only slightly changed as 
compared to the case of quantum dot in the absence of magnetic field, see Figs. |4] and 
[HJc). This is because Zeeman splitting does not provide a mechanism leading to a 
qualitative change of the noise - the noise and current are modified in a similar way, 
so the corresponding Fano factor is only weakly affected by magnetic field. 

5. Summary 

We have considered spin-polarized transport through a two-level quantum dot coupled 
to ferromagnetic leads in the sequential tunneling regime. The cases of symmetric and 
asymmetric coupling of the dot to external leads have been analyzed. For the most 
characteristic transport regions we have derived approximate analytical formulas for 
the current and Fano factor in the parallel and antiparallel configurations, as well 
as for the TMR effect. These formulas may be useful in the interpretation of future 
experimental data. 

Transport properties of quantum dots symmetrically coupled to the leads are 
shown to be strongly dependent on the bias and gate voltages. First of all, the bias 
and gate voltage dependence of TMR is then significantly different from that in the 
case of single-level quantum dots. The TMR can take now several well-defined values, 
which depend on the parameters of the system. Furthermore, we have shown that 
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the shot noise in the parallel configuration may become super-Poissonian in some 
transport regions. In the antiparallel configuration the noise is rather sub-Poissonian 
and the Fano factor becomes equal to unity when the leads become half-metallic. 

We have also shown that for a quantum dot coupled to one half-metallic and one 
nonmagnetic lead, the current is suppressed in certain bias regions. The suppression 
of the current results from the full occupation of the relevant dot spin state. A super- 
Poissonian shot noise has been found in these blockade regions. On the other hand, 
in the case when there is an asymmetry between the couplings of the dot levels to 
external leads, some blockade regions accompanied by NDC have been observed in 
both magnetic configurations of the system. We found three blockade regions with 
significantly different behavior of the TMR effect. In the regions where NDC occurs, 
the Fano factor is larger than unity and shot noise is super-Poissonian. 

Additionally, the effect of a finite Zeeman splitting of the dot levels was discussed 
in the case of quantum dot coupled symmetrically to the leads. We have shown 
that NDC occurs in both magnetic configurations, and for both bias polarizations. 
Furthermore, TMR exhibits peaks at voltages corresponding to consecutive Coulomb 
steps, and a high plateau in TMR appears in the linear response regime. Unlike the 
current and TMR, the Fano factor was found to be only slightly affected by an external 
magnetic field. 

The numerical results presented in this paper have been calculated for Ui = U2 = 
U' . However, the results for a general case, where the Coulomb integral between two 
electrons localized on different orbitals is different from the Coulomb integral for two 
electrons occupying the same orbital are qualitatively similar and the difference is 
rather of quantitative nature. More specifically, the lengths of different plateaus on 
the voltage scale may now be different, which is a consequence of the fact that the 
charging energy depends now on the distribution of electrons between the two orbitals 
of the quantum dot. Apart from this, the heights of the plateaus may also be different 
as compared to the case when all Coulomb integrals are equal. 
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Appendix 

In the following we describe the procedure to derive the analytical formulas presented 
in tableO As an example, we calculate the current in the parallel configuration for the 
transport region (3) shown in fig. 01 In this region there are five dot states taking part 
in transport, these are the empty state and all the four single-particle states, see table 
[1] As the occupations and contribution to the current coming from the other states 
are suppressed, it is sufficient to consider only the transitions (self-energies) between 
the above-mentioned five states. In the low temperature limit we approximate the 
Fermi-Dirac functions by step functions and get for the full self-energy matrix in the 
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basis defined by these five dot states 
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(p-i)r/2 
-(i+p)r/2 

{l-p)T/2 J 

r), which is due to normalization, one can 
calculate the occupation probabilities from Eq. ([3]). For plateau (3) and parallel 
configuration one finds that the probabilities are equal and given by 1/5, i.e. each of 
the five states is equally occupied. In a similar way one can determine the elements 
of self-energy matrix W'- 

/ -(i + p)r/2 (p-i)r/2 -(i+p)r/2 (p-i)r/2\ 

-(i + p)r/2 

= {p- i)r/2 

-(i+p)r/2 

\ (p-i)r/2 / 

It is now possible to calculate the current from Eq. which gives Ip ~ 2eT/'ah. 
The other analytical expressions can be determined in the same manner. 
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